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We study the ground state and the magnetization process of a spin-1/2 J1-J2 model with a plaquette structure 
by using various methods. For small inter-plaquette interaction, this model is expected to have a spin-gap and we 
computed the first- and the second excitation energies. If the gap of the lowest excitation closes, the correspond- 
ing particle condenses to form magnetic orders. By analyzing the quintet gap and magnetic interactions among 
the quintet excitations, we find a spin-nematic phase around J1/J2 —2 due to the strong frustration and the 
quantum effect. When high magnetic moment is applied, not the spin-1 excitations but the spin-2 ones soften 
and dictate the magnetization process. We apply a mean-field approximation to the effective Hamiltonian to find 
three different types of phases (a conventional BEC phase, "striped" supersolid phases and a 1/2-plateau). Un- 
like the BEC in spin-dimer systems, this BEC phase is not accompanied by transverse magnetization. Possible 
connection to the recently discovered spin-gap compound (CuCl)LaNb207 is discussed. 



I. INTRODUCTION 

Magnetic frustration have provided us with many intriguing 
topics e.g. the phenomena of order-by-disorder, the residual 
entropy at absolute zero temperature, disordered spin liquids, 
etoi. There are a variety of models which are known to exhibit 
the so-called frustration effects. Among them, the S* = 1/2 
J1-J2 model on a square lattice has been extensively investi- 
gated over the last two decades as one of the simplest models 
to study how frustration destroys magnetic orders and stabi- 
lizes paramagnetic phases. The model is defined by adding 
antiferromagnetic interactions on diagonal bonds to the or- 
dinary Heisenberg antiferromagnet on a square lattice (see 
FIG.IIJ: 



H ~ Ji ^ Si'Sj + J2 ^ Si-Sj 



(1) 



where the summations (n.n.) and (n.n.n.) are taken for the 
nearest-neighbor- and the second-neighbor (diagonal) pairs, 
respectively. In the classical (S / 00) limit, the ground state 
is readily obtained by computing Fourier transform J{k) of 
the exchange interactions and minimize it in the fc-space; 

• Ji > 0, J2 < Ji/2: the ground state has Neel antifer- 
romagnetic order (NAP). 

• J2 > \Ji\/'^- the ground state consists of two interpene- 
trating Neel-ordered square lattices. First quantum cor- 
rection fixes the relative angle between the two ordering 
directions and selects the so-called collinear antiferro- 
magnetic order (CAP). 

• otherwise: the ferromagnetic (PM) ground state is sta- 
bilized. 

Por J2 < 0, the next-nearest-neighbor (diagonal) interaction 
gives rise to no frustration and only the case with J2 > is 
non-trivial. The case Ji, J2 > has been extensively stud- 
ied in the context of spin-gap phases stabilized by the frus- 
trating interactions. Chandra and Doucot^ investigated the 
model in the large-S* limit and concluded that a non-magnetic 
(neither NAP nor CAP) phase appeared around the classi- 
cal phase boundary J2/J1 = 1/2. The most quantum case 



5 = 1/2 has been studied later both by numerical^ and by 
analytical methods^sS' (for other literatures, see, for instance, 
Refs. [tIS and references cited therein). By now it is fairly 
well established that we have spin gapped phase(s) in the win- 
dow 0.4 < J2/J1 < 0.6 although the nature of the spin-gap 
phase(s) is still in controversy. 

The case with Ji < 0, J2 > has been less investigated 
and recent analyses^ '" suggested that there is another non- 
magnetic (probably spin-nematic) phase around the classical 
boundary J2 / Ji ~ ^1/2 between CAP and PM. Prom an ex- 
perimental viewpoint, most compound a"i'^ found so far cor- 
respond to the ordered phase (CAP) of the J1-J2 model. 

Recently, Kageyama et al. reportedi^ a new two- 
dimensional Cu-based compound (CuCl)LaNb207. In 
this compound, two-dimensional sheets consisting of Cu^+ 
and CP are separated from each other by non-magnetic 
[LaNb207] layers and within each sheet the Cu^^ ions form 
a square lattice. The Cl~ ions are located at the center of pla- 
quettes and from a naive Goodenough-Kanamori argument the 
S = 1/2 J1-J2 model with Ji < and J2 > is suggested 
as the model Hamiltonian for (CuCl)LaNb207. 

What is remarkable with this compound is that inelastic 
neutron scattering experiment&i^ observed a finite spin gap 
2.3meV(=26.7K) above the spin-singlet ground state. Subse- 
quently, high-field magnetization measurementsii were car- 
ried out to show that magnetization monotonically increased 
between two critical fields Hd = 10. 3T and f/c2 = 30. IT. 
The data for (i) the Weiss temperature and (ii) the saturation 
field i7c2 in principle determine the coupling constants Ji and 
J2. Unfortunately, none of the solutions (Ji, J2) obtained in 
this way reproduced the spin-gap behavior—. Therefore, the 
usual J1-J2 model does not seem to work. 

The second intriguing point concerns the magnetization 
process. Prom the standard scenario'^, the onset of magne- 
tization at = i/ci in spin-gapped systems is understood as 
Bose-Einstein condensation (BEC, or superfluid onset, more 
precisely) of the lowest-lying triplet excitation (magnon) and 
the lower critical field i/ci at T = is given by the spin gap A 
as i7ci = A/{gfiB). This BEC scenario has been confirmed 
in various spin gap compound s 
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Recent specific-heat- and magnetization measurements^^ 
for (CuCl)LaNb207 exhibited behavior typical of spin-BEC 
transitions and suggested that the magnetization-onset tran- 
sition at may be described by BEC of a certain kind of 
magnetic excitations. However, we immediately find a seri- 
ous difficulty when we try to understand this within the stan- 
dard BEC scenario; the lower critical field H^i = 18. 4T 
expected from the observed spin gap A = 2.3meV at the 
zero field (where the experimental value g — 2.17 is used) 
in the standard scenario is much larger than the observed 
valuaii iJci = 10. 3T. One possible explanation for this dis- 
crepancy may be that a lower-lying triplet excitation which 
is responsible for the BEC was not observed in the neutron- 
scattering experiments because of selection rules. However, 
this seems unlikely since powder samples were used and usu- 
ally one can hardly expect a perfect extinction of a certain 
triplet excitation in such powder samples. Neither suscepti- 
bility measurementsii nor NMR data^° indicate such a hidden 
triplet excitation. 

An alternative and a more appealing scenario would be 
that the BEC occurs not in a single-particle channel but in 
a multi-particle channel. That is, what condenses to support 
a spin-superfluid is a bound state of magnon excitations. The 
possibility of multi-magnon condensation has been proposed 
theoretically2ii22 in the context of a kinetic quintet bound state 
in the Shastry-Sutherland model (see Refl23l and references 
cited therein). In fact, gapped quintet excitations which come 
down as the external field is increased were observed in the 
ESR experiments?"^ carried out for SrCu2(B03)2, whereas 
small Dzyaloshinskii-Moriya interactions hindered a quintet 
BEC from being observed in that compound (see also Refl25l). 

One of the simplest Ji-J2-like models which realize the 
above scenario and have a finite spin gap would be the S = 
1/2 J1-J2 model with a plaquette structure (see FIG. [U. A 
similar model (Ji, J2 > 0) has been investigated to develop 
a plaquette series expansion^. In this paper, we mainly focus 
on the region Ji < 0, J2 > where the quintet excitation is 
expected to play an important role in low-energy physics. 

The organization of the present paper is as follows. In sec- 
tion II, we briefly recapitulate the problem of a single pla- 
quette mainly to establish the notations. The coupling among 
plaquettes will be taken into account in section |III]by two dif- 
ferent methods; (i) a plaquette extension of the bond-operator 
mean-field theor}^ and (ii) a perturbation expansion with re- 
spect to the inter-plaquette couplings. We find gapped triplets 
and quintet for small enough inter-plaquette couplings in both 
methods. 

For larger values of inter-plaquette couplings, one of the 
gapped excitations softens and the form of the effective in- 
teractions among the soft excitations determines the resulting 
magnetic phases. By using the gaps obtained in the pertur- 
bation expansion, we determine the semi-quantitative phase 
diagram in section|lV](see FIG.fTOland FIG.[T2]|. 

The effect of high magnetic field will be considered in sec- 
tion |V] For high enough field compared with the spin gaps, 
we can approximate the low-energy sector by using only the 
singlet and the lowest excited state. For Ji < 0, we may 
expect that the quintet touches the singlet ground state first 



and a multi-particle BEC occurs. On general grounds, a 
single-particle (magnon) BEC phase is expected to have fi- 
nite transverse magnetization. Actually, in the BEC phase of 
TlCuCla, the transverse magnetization has been observed in 
the experiment!^. In the case of a multi-particle BEC, how- 
ever, the transverse magnetization does not appear To in- 
vestigate the magnetization process, we shall keep only the 
singlet and the quintet to derive a hardcore boson model as 
the effective Hamiltonian valid in high enough magnetic field. 
A mean-field approximation^ will be applied to the result- 
ing effective Hamiltonian to draw a full magnetization curve. 
Interesting phases (a 1/2 -plateau and supersolids) will be dis- 
cussed. According to the value of the parameters, we shall 
roughly classify the magnetization curve in FIG. [15] 

A summary of the main results and the discussion on the 
connection to the spin-gap compound (CuCl)LaNb207 will 
be given in sections |VI] and IVIII respectively. The equations 
omitted in the text will be summarized in the appendices. 

II. PLAQUETTE STRUCTURE 

We consider a spin- 1/2 J1-J2 model with a plaquette struc- 
ture where the interactions among spin-l/2s are explicitly 
tetramerized (see FIG. [U. The model is made up of four- 
spin units (plaquettes) and the four sites constituting a sin- 
gle plaquette are connected by the nearest-neighbor- ( Ji) and 
the second-neighbor ( J2) interactions as is shown in FIG. |2]l. 
The inter-plaquette interactions (both the nearest-neighbor- 
and the diagonal) which connect those units are multiplied by 
a distortion constant A (0 < A < 1). This parameter may be 
thought of as modeling the distortion of the underlying lattice 
in a simple way. In the case of A = 1, this model reduces to 
the homogeneous J1-J2 model, while when A = 0, the pla- 
quettes are decoupled from each other 





1 ( 

\ / 

/ s 

h — ^ 


I 




\ y 
s, ^ 

N/ 
/ >. 
/ S 

^ 


f V 


la 


\ y 
\ y 

/ \ 

f ^ 





FIG. 1: Two dimensional square lattice with a plaquette structure 
to be considered in this paper. Filled circles denote spin-l/2s con- 
nected by the usual exchange interactions. Thin lines (both solid and 
broken) imply that the interactions are multiplied by the distortion 
parameter A on these bonds. 



A. single plaquette 

Let us begin by analyzing a single isolated plaquette, which 
corresponds to the case A = 0. The eigenstates of a single 
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FIG. 2: Single plaquette. Dots represent spin l/2s, and the solid- and 
the dashed lines respectively represent Heisenberg interaction with 
the couplings Ji and J2 between S = 1/2 spins. 



plaquette can be easily obtained as follows. First we note that 
a plaquette Hamiltonian can be rewritten as 



H — Ji(Si-S2 + 82-83 + 83-84 
+ 72(81-83 + 82-84) 



84-81) 



4s' + ^(J2-Ji)(8^ + 8g)-^J2 



(2) 



where 8^ = 81 + S3, 8b = 82 + 84 and 8 = 8^ + 8^. There- 
fore, all the 2'^ eigenstates are classified by the three quan- 
tum numbers as |5'a, S*;,; S*). The eigenvalues E{Sa, Sb, S) 
are given by 



E{0,0,0) = , 
i;(l,l,0) = 2J2-2Ji, 

i;(i,o,i) = i;(o,i,i)- J2 

i;(l,l,l) = -Ji + 2J2, 

i;(i,i,2) = J1 + 2J2. 



(3a) 
(3b) 
(3c) 
(3d) 
(3e) 



Here a constant — ^J2 has been dropped just for simplicity. 
The energy of these states is shown in FIG. [3] For — 1 < 
J1/J2 < 0, the spin-singlet state |0, 0; 0) is the ground state, 
the triplets |1,0; 1), |0, 1; 1) are the first excited states, and 
the quintet |1,1;2) is the second excited state. For —2 < 
J1/J2 < —1, the singlet |0,0;0) is the ground state, quintet 
|1, 1; 2) is the first excited state, and triplets |1, 0; 1), |0, 1; 1) 
are the second excited state. 

The singlet |0, 0; 0) is written as 



|s) = |0,0;0) 

= ^(in)-UT))(in)-UT))- 

In what follows, the single-spin states in ket will be shown in 
the order of 1,3,2,4, i.e. |si, S3, S2, S4) in FIG. |2] For later 
convenience, we name the two triplets |1, 0; 1) and |0, 1; 1) as 
\pi) and \qi) (i ~ x,y,z) respectively. The explicit expres- 
sions of the two triplets are given as: 
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To label the quintet |1, 1; 2) states, we use the eigenvalues 
of S^, i.e. |1, 1; 5=2, S^) whose expressions are given ex- 
plicitly as: 



l,l;2,2) = |tTTT), 



(7a) 



i,i;2,i) = -{\ tT)(l Ti) + UT)) + (I Ti) + UT))I TT)} , 

(7b) 



1,1;2,0) 



V6 



{(I Ti) + I iT))(l n) + l iT)) + ITTii) + UiTT)}, 

(7c) 



1, 1; 2, -1) = -{| Ti) + I it)) + (I n) + I iT))l II)} , 

(7d) 

l,l;2,-2) = Uiii). (7e) 
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FIG. 3: The energy of the triplets |1, 0; 1), |0, 1; 1) and the quintet 
|1, 1; 2). We take the units of energy as J2, and the energy is plotted 
as a function of Ji / J2 



III. EFFECT OF INTER-PLAQUETTE INTERACTION 

For A = and J1/J2 > —2, all plaquettes are in the singlet 
state |0, 0; 0). Finite inter-plaquette interactions A induce var- 
ious tunneling processes among plaquettes to change both the 
ground state and the excitations over it. For finite A, we calcu- 
late the excitation energy by two different approaches. One is 
the bond-operator mean-field theory (MFT)2^^, which gives 
the excitation energy of the triplets \pi), \qi). Another is the 
second-order perturbation theory in A, and it gives the energy 
of the quintet |1, 1; 2) as well as that of \pi) and \qi) . For 
sufficiently small A, both approximations yield finite energy 
gaps for these excitations and when one of these gaps closes, 
the corresponding (bosonic) excitation condenses to form a 
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magnetically ordered state. The energy of triplet excitations 
can be observed by inelastic neutron scattering experiments. 
Both approximations may not be reliable for large A and small 



A. bond-operator MFT 

Let us begin with the bond-operator MFT^^i^. For — 2 < 
•h/ J2 < and A = 0, |0, 0; 0) is the ground state and the 
degenerate triplets |1, 0; 1), |0, 1; 1) are the first- or the sec- 
ond excited state (see FIG.|3]l. Therefore, we may truncate the 
Hilbert space and consider a subspace spanned by the singlet 
|0,0;0) and the triplets \pi), \qi). This approximation is re- 
liable to estimate the excitation energy of the triplets, unless 
1 1, 1; 2) condenses. In this subspace, nonzero matrix elements 
of Si, 2, 3,4 is 





2^al3 , 


{Pa 


S?\Py) = 


i 

2^at3j 7 


(8a) 






{la 
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(8b) 


s\SS\Pf3) = 


~2'^"'3 ' 


{Pa 


5f b7> = 


i 

2^al3f 1 


(8c) 






{qa 




i 

2^af3f 7 


(8d) 



where a, l3—x,y, z. Using boson operators s,pc, (7^(0; = 
X, y, z) satisfying the standard commutation relations, 

[SjSl'] = 1, [Pa,p'^p] = [(la,q\] = (5q/3, [s,Pa] = 0, 

etc, the local spin operator Si 2,3.4 can be written as 



1 i 

-(s^p^ + spt^) _ g^^^ptp^ ^ (9b) 



2' 2 

'^3-^2^""" ' 2^ 



(9a) 



7 I 



(9c) 
(9d) 



where the summation over repeated indices is implied. Since 
the restriction that each plaquette has exactly one particle 
leads to the local constraint s^s + J^aiPaPa + qt^a) = 1, 
we introduce the Lagrange multiplier /i; and add a term 



^J■^ 



Y.^pI,zPo 



qhqo 



(10) 



to each plaquette Hamiltonian. We may assume that //^ for 
each plaquette takes the same value ji for all plaquettes be- 
cause of the translation invariance. 

Next, we replace s by its expectation value (s) = s, since 
the s boson condenses in the ground state. Moreover, since 
the triplet is dilute when the energy gap is positive, we may 
ignore the terms consisting of three or four triplet operators. 
In this way, we obtain the mean-field Hamiltonian i^bo con- 
sists only of bilinear terms in p and q. The mean-field param- 
eters {fi, s) are determined by requiring the expectation values 



of the derivatives of _Hbo with respect to the mean-field (MF) 
ground state vanish: 



(f^) =0, (^) ^0, (11) 

\ Ofl / MF \ OS / MF 

or equivalently by finding the extrema of the mean-field 
ground-state energy E'q ^ : 



/^H^ 



dfjL 







(12) 



In particular, Eq\ must be minimum for s. 

In this approximation, the inter-plaquette interactions asso- 
ciated with the site n reads 

{,Hx)n = ^l(S2-Sa + Ss'St) + J2(S3-Sa + S2-S5) 
= '^'s'^{{qa + qi)n[Pa + P^a)n+x + {Pa + p\)n{qa + qi)n+x} 
- '^'s'^{{Pa + Pa)n{Pa+Pa)n+x + {qa+qi)n{qa + qi)n+x), 

(13a) 

{Hy)n = >/l(S3-Srf + S4-Sc) + J2(S3-Sc + S4-Sd) 



S^{{q» + qa)n{Pa + Pa)n+y+iPa + Pa)n{qa + qi)n+y} 



J2-2 



s'-'{{Pc, + Pl)n{Pc. + p'^a)n+y + {qa + qi)n{qa + qi)n+y}, 

(13b) 
(13c) 



{Hx+y)n — J2^3'^g 

= -^S^(Pa +Pi)n{Pa+pi)7i+x+y, 



{.Hx-y)n — J2S2-S/ 

= -^s'^iqa + qi)n{qa + qi)n+x-y 



(13d) 



where the site labels 1, . . . , 4 and a, . . . ,g are defined in 
FIG.H 
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1^ Q- 
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FIG. 4: Inter-plaquette interactions associated with the plaquette n 
(shown by a thick line). 



Summing up all four interactions and doing Fourier trans- 
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formation, the total Hamiltonian i?bo reads 

k 



4 

J2-2 

4 " 



s /-(k)(gk 9k + 9k9k + 9k9-k + 9k 9-k) 

2 „ 



-M^'+pi.>k + 9l"9£-l)}, (14a) 
where we have defined 

/±(k) = cosfcj; + COS fey + cos{kx ± ky) . (15) 

If we introduce a vector Vk = (pk , 9k , pLk i 9-k)^' '■^^^ 
MF Hamiltonian iJbo can be written compactly as 

= 5] vt"A(k)v£ - iVp {3(J2 - m) + - 1)} , 

k 

(16) 

where A'p denotes the total number of plaquettes and the ker- 
nel A(k) is given as 



A{k) = 



/ abo(k) &bo(k) Cbo(k) 4o(k) \ 

6bo(k) ebo(k) &bo(k) 4o(k) | 

Cbo(k) 6bo(k) abo(k) 6bo(k) 

V6bo(k) 4o(k) 6bo(k) ebo(k)/ 



(17a) 



abo(k) = 


'h J2-2 f n \ 


(18a) 


&bo(k) = 


— ^s^{coskx — cosfcy) , 


(18b) 


Cbo(k) = 


-f -f+w ■ 


(18c) 


4o(k) = 




(18d) 


ebo(k) = 




(18e) 



Using a 4 X 4 real matrix Lk (see Appendix|A]for the detail), 
we can diagonalize ^(k) by the Bogoliubov transformation: 



(19) 



ikVk = Vk , 

'flip t' 
Vk = 1^ Pk ' 9k ' P-k ' 9-k 

As is shown in AppendixlAl i/bo then reduces to 

i^bo = J2 {^i(kW>L" + u;2{k)q'^"qL"} + E^i , (20) 



where the mean-field ground state energy is given as 
3 



-(c^i(k) + i02ik)) - 3(J2 - 1) 



(wi(k) , cJ2(k)) = (cj(+._)(k) , W(+.+)(k)) , 



W(±,±) ^ 

+ 4(a 



±2 

-d)(c 



(21a) 
(21b) 



(21c) 



In eq. ( |21c| ). the order of signs ± coincides on both sides. 
Since 0-11^2 > 0, condensation of the triplets p and q occurs 
when the equality holds at some k. Otherwise, there is no 
condensation, and (pjf ) = WC) ~ 0- Therefore, there exist 
rotational symmetry and no magnetic order In this case, cji_2 
are the excitation energy of triplets. 

We looked for the solutions (^,s) to the set of equa- 
tions (fT2l i numerically. For example, we found (/i,s) = 
(-0.09,0.96) for the set of parameters A = 0.3, J1/J2 = 
—0.8. The dispersion relation of the excitation energy 
u!i(k) =W(+__)(k) is shown in FIG. |5] 




no. 5: The dispersion relation of the excitation energy of the triplet 
p' in ( 120b . which has the lower energy of the two triplets, for the 
parameters A — 0.3, Ji — —0.8, J2 = 1 



If the excitation becomes soft w = at some k, the sys- 
tem is in a magnetically ordered phase. From the known 
results^iiS, we expect that ordered phase appears for A suffi- 
ciently close to 1 . To determine the phase boundary between 
the paramagnetic phase and magnetically ordered ones, we 
searched the (A, Ji / J2) plane for the points where the mean- 
field gap vanishes. Unfortunately we found that the gap did 
not close in the relevant parameter region 0<A<1, — 2< 
•h/ J2 < 0, and that the disordered singlet phase persisted; the 
gap vanished only for larger A(> 1). This unacceptable result 
may be attributed to the fact that the bond-operator mean-field 
theory probably overestimates the stability of the plaquette 
phase. 
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B. Second order perturbation 

In this section, we compute the energy gap of triplets 
|1,0; 1), |0, 1; 1) and the quintet |1, 1; 2) by the second or- 
der perturbation theory in the distortion parameter A. The 
naive expansion in A is ill-behaved in the vicinity of the point 
Ji/ J2 = —2 and we have to use another perturbation scheme 
for that region. 



1. The excitation energies of triplets 

Let us consider the states where there exists only one triplet 
and all the other plaquettes are in the singlet |0, 0; 0) state. If 
the coupling constant of inter-plaquette interaction A=0, these 
states are iVp-fold degenerate, where iVp is the number of pla- 
quettes. For finite A, the second order perturbation induces 
hopping of the triplet to nearest or next nearest neighbors and 
lifts the degeneracy. 

Rotational symmetry forbids the hopping which changes 
the spin label i{~ x,y,z) or the magnetic quantum num- 
ber. On the other hands, the transitions between two different 
triplets Pi and qi of the same label i occur For example, the 
hopping amplitude of Pi{qi) to the nearest-neighbor plaquette 
is given by 



given by eq.( ID4l ). The dispersion relation of the lower branch 
is shown in FIG.|6] 




FIG. 6: The dispersion relation of the excitation energy (k) of 
triplet at A = 0.3, Ji = -0.8, J2 = 1. 



(22) 



The degeneracy is partially resolved by the hopping of pi{qi). 
The transition between pi and q, will be considered later. In 
the second-order perturbation, the processes that the triplet re- 
turns to the original site is also allowed. Including this effect, 
the energy change of pi(k)-particle is given by a(k) in (IDlb . 
Similarly, that of qi (k) is given by 5(k) in ( ID2I ). 

Next, we consider the transition between Pi{k.) and gi(k). 
The transition amplitude is given by 



The lower branch takes its minimum at the F-point 
k = 0, and E^ (k — 0) gives spin gap A^. The second or- 
der expression of At is given in ( ID5I ). The expression tells 
us that At has a pole at Ji = — 2J2 and that the standard 
perturbation breaks down near the pole. To remedy this, we 
introduce another perturbation parameter 5 = Ji — (— 2J2) 
and carry out a double expansion in both A and S. Then, we 
obtain the energy gap in a modified method E^^^^^{0) given 
in eq.( ID6l ). This improved energy gap is expressed to give a 
better approximation around Ji = — 2 J2 



c(k)^ 



2^+4J2^ 



(cos — cos ky) 



(23) 



Therefore, for each i = {x, y, z), eigenstates t+,t- satisfies 



2. excitation energy of quintet 



a(k) c(k) 
c(k) 6(k) 




(24) 



The expressions of a(k), fe(k), c(k) are given in Appendix Id] 
After this procedure, the degeneracy with respect both to the 
position and to the species pi and qi is resolved. There also 
exists an energy shift in the ground state. Taking all these into 
account, we obtain the energy of the triplets: 



^ ^a(k)+5(k)± V(«(k)-6(k))2 



- 4c(k)' 



(25) 



where Ai?s denotes the energy shift of the bare ground state 
where all plaquettes are occupied by the singlet |0, 0; 0) and is 



Next, we consider states containing only one quintet in a 
background of the singlet plaquettes. As before, the degen- 
eracy with respect to the position of the quintet plaquette is 
resolved by hopping. Up to the second order in A, the hop- 
ping to nearest neighbor is given by 



v2 ^1 



J? 



8J1 



(26) 



and the hopping to next nearest neighbor does not occur. Tak- 
ing into account the processes that the quintet returns to the 
original site and the energy shift of the ground state, the exci- 
tation energy of quintet is given by Eq{k) in (|D7| i. The dis- 
persion relation is shown in FIG. |7] 



FIG. 7: The dispersion relation of the excitation energy -Bg(k) of 
quintet at A = 0.3, Ji = -0.8, J-z = 1. 



J1/J2 

FIG. 8: The value of A at which the energy gaps At and Aq close. 
The energy gap of the triplets At is given in l IDSt . which is not re- 
liable near J1/J2 = —2 because of the pole there, and that of the 
quintet Aq is in l l27b . which is not reliable near Ji / J2 = 0. 



Since the quintet dispersion i?g(k) takes its minimum at 
k = 0, the quintet gap is given by 



A, = i;,(k = o) 

= Ji + 2 J2 - 



A^(6Jf-23JiJ2 + 10Jf) (27) 
8Ji(2Ji - J2) • 



We note that there is the pole at Ji = and the approximation 
becomes poor for Ji w 0. 



IV. GROUND STATE PHASES 

If the inter-plaquette coupling A is increased, one of the en- 
ergy gaps of the triplets (( ID5I ) and ( ID6b ) and the quintet ( ID7b 
becomes at a certain critical value of A. When it happens, the 
corresponding particle condenses and a phase transition oc- 
curs from the gapped spin-singlet phase to superfluid phases 
with magnetic long-range order. Therefore, we can classify 
the phases according to what kind of particles condense and 
what kind of magnetic orders is stabilized by a given set of 
interactions among them. In FIGs.[8]and|9l we plot the value 
of A at which the smallest energy gap becomes 0. If we as- 
sume that no further condensation occurs in the other kinds of 
particles once the triplets or the quintet condenses, the phase 
diagram FIG.fTOlis obtained. When we mapped out the phase 
diagram FIG.[TOl we have used two different expressions ( |D6b 
and dPSb for the energy gap of the lowest triplet in the vicinity 
of J1/J2 = —2 and away from it (J1/J2 ~ 0), respectively. 
We have also neglected the quintet around Ji / J2 = since 
the collapse of the quintet gap there (see FIG. [8]) can be at- 
tributed to the existence of a pole and is just an artifact of the 
perturbative approximation. Note that the phase boundary be- 
tween the two regions covered by eq. dPSb and eq.( ID6l ) is only 
schematic. 

Now let us discuss the nature of the ordered phases appear- 
ing after the condensation. In the region shown as "CAF" 




-2 -1 .95-1 .9-1 .85-1 .8-1 .75-1 .7-1 .65-1 .6-1 .55-1 .5 

J1/J2 

FIG. 9: Plot of the value of A when the energy gaps are equal to 0. 
We use At, mod (eq.( |D6b ) for the triplets, which is reliable even in the 
vicinity of J1/J2 = —2. The curve for the quintet is the same as in 
FIG.[8] We only plot the region -2 < J1/J2 < -1.5. 



(highlighted in red) in FIG. [TOl condensation occurs to the 
singlet and the triplets. Then, we may expect: 

\{s)\ ^ 0, |(p(k = 0))|2 ^ 0, |(q(k = 0))|V , (28) 

which, combined with ( |9dl i. implies 



(51) = -(83)^0, 

(52) = -(S4) ^ , 



(29a) 
(29b) 



provided that e"'37(pt)(p^) ^ Q a.nA e'^P^ {q^) {q-,) = 0. Note 
that all the plaquettes are in the same state, since the energy 
of the triplet takes its minimum at the F-point k = (see 
FIG. m. When the combination (p + q) of the two bosons 
condenses, the relation (Si) = (S2) holds and the ground 
state has the transversely aligned (i.e. (0, it)) collinear anti- 
ferromagnetic order In the case where (p — q) condenses, on 
the other hand, we have (Si) = — (S2) instead and the system 
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FIG. 10: The schematic phase diagram of the ground state deter- 
mined by the particle whose excitation gap closes first. In the green 
region, the quintet and the singlet condense, in the red do the triplet 
and the singlet, and in the blue does the singlet. In the region marked 
by blue, the energy gap exists. The phase shown by red may be con- 
sidered as coUinear antiferromagnetic (CAF) state. The nature of the 
green phase, where the quintet condensation occurs, is closely inves- 
tigated by using an effective Hamiltonian H^^ (eq.l|32|l). 



is in the collinear antiferromagnetic ground state in the longi- 
tudinal ((tt, 0)) direction. This is consistent with the known 
results^. 

Now we move on to a more interesting case. In the green 
region in FIG.fTOl frustration is strong ( J2 ~ — 2 Ji) and non- 
trivial order may be expected. In fact. Shannon et al^^ ana- 
lyzed the uniform (A = 1) model by numerical exact diago- 
nalizations up to clusters of 36 spins and found a spin-nematic 
phase with d-wave (or Bi) symmetry for —2.5 < J1/J2 < 
— 1.43 ~ —1.67. In the state with the nematic order, the ex- 
pectation value of the rank-1 tensor vanishes (Si) — 0, while 
we have a finite expectation value of the following traceless 
rank-2 tensor: 

Ql' ^ ' - - ^S'^' , (30) 

where o;, /3 = (x, y, z) and i.j label the lattice sites. 

As is shown in FIG. [TOl the singlet and the quintet con- 
dense in the region of interest. This is analogous to the spinor 
Bose-Einstein condensation of spin-2 particles (here particles 
are defined not on the lattice sites but on the plaquettes). We 
consider a single plaquette (see FIG.|2| and, as before, denote 
the singlet and the quintet respectively by | s) and 1 1 , 1 ; 2 , 6*2) . 
To investigate what kind of magnetic order is stabilized in the 



condensate, let us introduce the following mean-field ansatz 
for the ground state: 

im,{V'r}) = 

(g) JcOS0r|s)r+sin0r ^ MS z)\l S A , 

replaq I S^ = -2 ) 

(31) 

where the product is over all plaquettes and the complex num- 
bers tpr{Sz) satisfy J^s H'r{Sz)\'^ = 1- Then, since therank- 
1 tensor can not give rise to transitions between the spin-0 
states and the spin-2 ones by the Wigner-Eckart theoremi-1, 
we have (5|Sf|s) = {s\Qf\s) - 0, (sj^f |1, 1; 2, 5z) = 
and consequently ({6ir},{V^r} I S, I {6ir},{i/;r}) = 0. If 
we introduce the cyclic operator C which translates the 
state as 1 ^ 2 ^ 3 ^ 4 ^ 1, we obtain 
C\s) ^ -\s), C|l,l;2,S'z) = |l,l;2,S'z) from © 
and Therefore, the spin-nematic tensor Q"^f^ de- 

fined on the bond satisfies (s|(5"|'|l, 1; 2, S'z) = 

(s|CtCQ?2^CtC|l,l;2,5z) = -{s\Q^^\l,l;2^ Sz) = 
{s\Q'^!^\l,l;2,Sz) = -~{s\Q2^\l,l-2,Sz). This implies that 
the spinor condensate \{0r}, {^r}) of our quintet boson has 
the same (d-wave) symmetry as the spin-nematic state dis- 
cussed in Ref|Tci. 

However, this is not the end of the story. Since the 
local spin operator with S > 1 assumes several dif- 
ferent states (e.g. polarized, nematic, etc.) and it is 

not obvious if {{Or},{i'r}\Qtf\{Or},{M) ^ OI UOt 

for our J1-J2 model. To determine the actual value of 
{{Or} , {'4'r}\Q'if \{(^r} , {i/Jr}) , wc need the explicit mean- 
field solution for a given set of (Ji, J2,A). Since we are 
considering the situation where the gap between the singlet 
ground state and the quintet excitation is vanishingly small, it 
would be legitimate to keep only the singlet |0, 0, 0) and the 
quintet for each plaquette to write down the effective Hamil- 
tonian. 

The form of the effective Hamiltonian is determined by us- 
ing the second-order perturbation theory and it contains the 
kinetic part describing the hopping of the quintet particles 
and the magnetic part which concerns the interactions among 
them. Since within a mean-field treatment the spinor part 
i-'riSz) is determined by the magnetic interactions, it suffices 
to consider only the magnetic part of the effective Hamilto- 
nian: 



{i".j".k") (32) 
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where S"? denotes the 5 = 2 spin operator, and the symbols 
and {i',j') mean the nearest-neighbor- and the next- 
nearest-neighbor pairs, respectively. For different types of 
three-plaquette clusters k"), we assign different three- 

body (i.e. three-plaquette) interactions L)^^!^ 2 = 1j • ■ • 1 6) 
in ( l32b . The correspondence between six types of clusters and 
the strength of the three-plaquette interaction L^^j^ 2 is shown 

in FIG, nn The full expressions of Jqui,2, ^i^qui,2 and L^"^ 2 
are given in Appendix. |Dl Note that our effective Hamiltonian 
in its full form contains the kinetic term and charge interac- 
tions as well as magnetic ones i/qu- In this sense, our effec- 
tive model is a generalization of the Bose Hubbard Hamilto- 
nian for F = 2 cold atoms in optical Iattices22ii2i21 and the 
determination of the full phase diagram and the identification 
of various phases found in systems of cold atoms in our mag- 
netic system would be interesting in its own right. 

~ L 

1 2 3 

4 5 6 

FIG. 11: Clusters involved in the 3-points interaction in ( i32t . The 
plaquette corresponding to i" is always located on the center of the 
clusters. We identify all clusters obtained from a given one by rota- 
tion by it/2, tt, 3n/2 and reflection. 



We investigate this Hamiltonian by means of a mean-field 
theory by assuming an r-independent uniform {0, V'}, for sim- 
plicity. Since the parametrization of the spin-2 states is cum- 
bersome, we adopt the method used by Bacry^- and Barnett et 
alJ^. First we note that arbitrary (normalized) spin-S" states 
are parametrized by a set of 25 unit vectors except for ob- 
vious gauge redundancy. Using rotational symmetry, we can 
further reduce the number of free parameters needed to ex- 
press arbitrary spin-2 states to 2x4—3=5 (see Appendix. IbTi. 
We numerically minimized the mean-field energy with respect 
to these five parameters. The result is shown in FIG. [12] 

At A = 1, the system is in the ferromagnetic state for 
•h/ J2 < —2.33 and is in the spin-nematic state for —2.33 < 
Ji/J2{< —1.91). This result slightly differs from the numer- 
ical results2i£ -2.5 < J1/J2 < -1.43 1.67. However, 

this is not surprising since our results are based on a mean- 
field treatment of the magnetic Hamiltonian Hq^ obtained by 
perturbation expansion in A. Our result may be improved by 
taking the number of sublattice larger, since Jqu2 > and 
there are various 3-site interactions Lq^i and Lqu2- 



-233 -1.91 




I I I 

-2.6 -2 -i.5 

J1/J2 

FIG. 12: The schematic phase diagram obtained in a similar man- 
ner to in FIG. [To] We zoom up the region around J1/J2 — —2 in 
FIG.[TO] In the two regions on the left (green and yellow), the quintet 
and the singlet condense and we determined the resulting magnetic 
orders by a mean-field approximation to the magnetic Hamiltonian 
Hqu. In the green region, the quintet and the singlet condense, and 
the spin-nematic phase appears. In the red, on the other hand, con- 
ventional ferromagnetic order is stabilized. The red and the blue 
region are the same as FIG. 1 101 



V. MAGNETIZATION PROCESS 



Having mapped out the phase diagram in the absence of 
external magnetic field, we consider next the magnetiza- 
tion process of the plaquette model by mapping the orig- 
inal model onto a hardcore boson model or an equivalent 
5 = 1/2 pseudo-spin model. Tachiki and Yamada^ applied 
this method to obtain the magnetization curve of the spin- 
dimer model, which consists of pairs of 5 = 1/2 spins. The 
coupling to the external magnetic field is incorporated into the 
Hamiltonian by adding the Zeeman term g/ish- J2i ^i- For 
convenience, we set g/ie = 1 and assume that h is pointing 
the z-direction: h = (0, 0, h). 

Although the original treatment in Ref. [s^ is for a coupled 
dimer systems, we can readily generalize the method to our 
plaquette system as follows. We denote the plaquette states 
by |5q, 5fc; 5, 5^), where 5, 5q, St are defined in ((2|. From 
([3]i, the energies of a single plaquette satisfy 



S(l,l,2)<2i;(l,0,l), 2i;(0,l,l), (33) 



for — 2 < Ji/ J2 < 0. As is shown in FIG. [13] with increasing 
the magnetic field, the quintet level |1, 1; 2, — 2) comes down 
to |0, 0; 0) faster than the lowest triplet levels |1, 0; 1 — 1) and 
|0,1;1,-1). 
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E 
A 

Ji+2J. . 

\ ^ 1.0;l.-i> 
|l,l;2,-2> 

FIG. 13: The energy of eigenstates of a single plaquette as a function 
of magnetic field h. 



spanned by |0,0;0) and |1,1;2,— 2) probably breaks down 
for ft, « where all the components (S^ = — 2, . . . , 2) of the 
quintet come into play. Also the validity of the approximation 
may be questionable for sufficiently large A where the singlet- 
triplet gap may be much smaller than the singlet-quintet gap, 
since the triplet states |1, 0; 1) and |0, 1; 1) are important there 
(seeFlG.Uni. 

If we simply project the original S" = 1/2 Hamiltonian 
to the restricted subspace as in ([8]), no spin-flipping term (or, 
hopping term, in terms of hardcore bosons) appears. This is 
because the projection is equivalent to the ordinary first-order 
perturbation theory and no transition between the singlet and 
the quintet occurs in the first-order processes. Therefore, we 
need take into account the second-order processes to obtain 
the meaningful effective Hamiltonian. The amplitude that a 
quintet state (spin 'down') |1, 1; 2, —2) hops to the adjacent 
plaquette is given by 



Therefore, in order to describe the low-energy physics in 
the presence of strong magnetic field (h Ji/2 + J2), 
we may keep only the two lowest-lying states |0, 0;0) and 
1 1,1; 2,— 2) for each plaquette and restrict ourselves to the 
subspace spanned by them. In what follows, we regard the 
singlet |0, 0; 0) and the quintet |1, 1; 2, —2) respectively as the 
up- and the down state of a pseudo spin- 1/2. That is, 

|0,0;0)= Q , |l,l;2,-2)= (^^ . (34) 

Then, the resulting effective Hamiltonian is written in terms of 
the Pauli matrices {S = 1/2 spins) defined on each strongly- 
coupled plaquette. 

Note that the approximation to treat only the subspace 



72 I 72 

t^X'^^l±^. (35) 

The hopping to the next nearest-neighbor does not occur at 
this order of approximation. The energy gap between the state 
where there exists only one static 'down' spin (|1, 1; 2, — 2)) 
in a background of the 'up' spins (singlet |0, 0; 0) plaquettes) 
and the one where all plaquettes are 'up' is given by —fi in 
( ID13ab . The interaction between the two adjacent 'up' spins 
(|1,1;2,— 2))is given by Joffi in ( ID13bb and that between the 
next-nearest-neighbor pair is given by Jcfr2 in ( ID13cb . We 
note that this approximation becomes poor near the pole of 
•/cff i,2 and t at J1/J2 ~ 0. On top of them, we have several 
three- 'site' processes and putting them all together, we obtain 
the effective Hamiltonian: 



i?eff = ( Ji + 2J2 - - 2h) <y~^t + E {^('^t^J + '^r^) + '/effi(fT-a+)(a-t7+)} 

+ -^-ff^ (^^^ ) 4 ) + J2 ) 4' ) (-^fe" ' ) ' 

(i" ,j" ,k") 

I 



where crs denote the Pauli matrices and cr^ = ^(cr^ + 
1(7^), cr^ = ^(cr^ — icr^). The symbols and 
mean that the summation is taken over the nearest-neighbor- 
and the next-nearest-neighbor plaquettes, respectively. As in 
section HV] there are six types of Loff for different bond con- 
figurations k") (see FlG.[n]l. We label the different 

three-plaquette interactions by L^g'' (n = 1 . . . 6) and the cor- 
responding bond-configurations are shown in FIG. [TT] The 
concrete expressions of Jeff and igff are given in Appendix. ID] 
We note that the transverse components ct+ and can be 
translated to the creation- a and the annihilation operator 
of a hardcore boson, respectively. 



We analyze the Hamiltonian (|36] | within a mean-field ap- 
proximation. Since Jcffi.2, which have the first order contri- 
butions in A, are dominant for small A, we may assume two 
different two-sublattice structures: (i) "checkerboard" and (ii) 
"stripe" shown in FlG.[T4lin the calculation . 
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we can rewrite ( 1361 1 in terms of (T*(i = x,y,z). Since we 
are interested in the ground state energy at T = 0, we can 
simply replace operators in (|36] | by their expectation values 
on each site, e.g. ^ E(zj> (c^^) (o""") for the 

"checkerboard" case. For convenience, we introduce the fol- 
lowing two-component vector: 



FIG. 14: Two-sublattice structures assumed in the calculation: (i) 
striped- (left) and (ii) checkerboard (right) case. Circles (whether 
filled or open) denote the strongly-coupled plaquettes shown by thick 
lines in F1G.|I] 



(38) 



By using the relations 

a-a+ = ^{l-a') , (37a) 
ara+ + a+aj = i^a^ + af aj) , (37b) 



Since there is rotational symmetry in the x-y plane, the mean- 
field energy is parametrized by (cr^), (cr^'), r = |r|, t' = \t\ 
and the angle </> between r and t' . The Hamiltonian (l36l l 
reduces to 



+ /3irT C0S(/) + /?2^ +^1- + 12- t: 



r 



(39) 



where TVp denotes the total number of plaquettes and a, /3, 7, ^ 
are given in Appendix |D] both for the case of "checkerboard" 
and for the "striped" case. Correspondingly, the total magne- 
tization is given simply as 



M 



— E 

2iV„ ^ 



(40) 



In both cases, /?i 2 < and E^s is minimized for = 0. 
Since any spin- 1/2 states satisfy the following relation among 
the expectation values (see eq. ( IC5I )) 



1 



(41) 



the transverse magnetization r can be expressed in terms of 
the longitudinal one (cr^). Hence, there remain two variational 
parameters (cr^) and (u^') in £^eff- From the definition ( l34l i. 
the expectation values (cr^) = 1 and (tr^) = — 1 respectively 
correspond to the singlet state and the fully polarized (or, sat- 
urated) state. 

The critical field h ~ H^i which marks the onset of mag- 
netization is given by (9£'cff/9(cr^))(o-^)=i = after substi- 
tuting (<T"') = (cr^),i.e. 



2Hci = Ji + 2 J2 - M + 4i . 
The right-hand side is exactly the same as ( |27] |. 



(42) 



Once spin-gap closes at h = H^i, the quintet particle 
|1, 1; 2, -2) condenses, i.e. {a'') 7^ 1, r ^ 0. If (cr^) ^ ±1, 
r ^ and there exists a finite expectation value of (cr^). 
In the hardcore boson language discussed below (|36] |, 
can be viewed as the boson annihilation operator a and its 
finite expectation value (a) ^ implies that Bose-Einstein 
condensation of the quintet particle occurs. In particular, if 
I ('^^) I 7^ I ('^^') I ^iid T 7^ r' in BBC phase, the state is in the 
so-called "supersolid" phaseii. For convenience, we shall call 



the BEC phase satisfying (cr^ 



(cr^') a normal BEC. 



It should be noted that even when r 7^ 0, the transverse 
magnetization (S*^) vanishes unlike the BEC in the spin- 
dimer modeU^. In fact, since the creation operator at of the 
quintet particle can be written in terms of the original spin 
operators as 



= ^jqb:-qii)+^qi', 

Bi ^12 ^ W23 T" "^34 ~ ^41 J 



(43) 



the existence of the condensate (cr+) 7^ (or t 7^ 0) implies 
that we have a finite expectation value of the following spin- 
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nematic operator: 



{Qi:-Ql\)±^Q 



(44) 



The form (03]) of the quintet creation operator suggests that 
we should think of the plaquette quintet 1 1, 1; 2, 2) as a tightly- 
bound magnon pair (or magnon molecule). 

The critical field i/c2 where the saturation occurs is given 

by {dE^s/d{(T''))i^„.)^_i = after substituting (tr^') = 
(a^),i.e. 



2Hc2 — Ji + 2J2 



4<- 



9r(l) or(3) o (4) 



4Jeffl + 4 Jeff2 
(6) 



2L, 



(45) 




eff 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 

FIG. 16: Magnetization curve for various values of the distortion pa- 
rameter A. The frustration parameters is fixed to Ji/ J2 = —1.4. The 
colors of the curves correspond to those used in FIG. [T5] (except for 
yellow). The blue curve has supersolid phase around the 1/2-plateau 
and the phase transition between the normal BEC and the supersolid 
phase is of second-order. All curves in BEC and supersolid phase is 
convex down because of the 3-point interaction 7 in l |39l ). 



We minimized E^s numerically and we found that the energy 
in the "stripe" case was always equal to or smaller than that 
in the "checkerboard" case. We show various types of mag- 
netization curves obtained in this way in FIG. [16] In FIG. [15] 
we also classified the parameter regions (in the {Ji/J2,X)- 
plane) according to the qualitative behavior of the magne- 
tization curve. There appears (i) the normal BEC phase, 
(ii) the "striped" supersolid phase and (iii) the "striped" 1/2- 
plateau. At the 1/2-plateau, the pseudo-spins a are ordered in 
a collinear manner (cr^) = 1 and (cr^') = —1 (see FlG.[T4b. 



I 



A 




FIG. 15: Schematic classification of the magnetization curve, (i) 
In the green region, the curve is smooth and the system is always 
in the normal BEC phase, (ii) In the red region, the magnetization 
curve has a 1/2-plateau. Except at the plateau, the system is in the 
normal BEC phase, (iii) The region where we have additional su- 
persolid phases around the 1/2-plateau is highlighted in blue, (iv) In 
the region colored by yellow magnetization jumps to saturation and 
the magnetization process is step-like. The concrete expression of 
curves is shown in FIG.[T6] The phase boundary is only schematic. 



The magnetization curve in the BEC and the supersolid 
phase is convex down because of 3-point interaction 7 in ( |39l ) 
which breaks the particle-hole symmetry. The "striped" su- 
persolid phase always appears around the 1/2-plateau and the 
width of the supersolid phase appearing on the left of the 1/2 
plateau is broader than that on the right due to the convex 
down character The equivalent Hamiltonian ( |36l ) without the 
3-point interactions has been investigated by using the mean 
field theory^ and Monte-Carlo simulation s^-' - " '. They found 
that the "striped" supersolid phase around the 1/2-plateau is 
stable^''. Therefore, our result that the supersolid phase exists 
may be correct beyond the mean-field approximation, since 
the 3-point interaction in ( |36] l is weak. There are other mod- 
els accompanied by the supersolid phase, e.g. spin dimer XXZ 
model^i, spin- 1/2 XXZ model on the triangular lattice^, etc. 



VI. COMPARISON WITH THE EXPERIMENTAL DATA 
OF (CuCl)LaNb207 

In this section, we compare our results with the experimen- 
tal data obtained for (CuCl)LaNb207. Since we have three 
parameters Ji,J2 and A, three experimental inputs in princi- 
ple determine the set of coupling constants. Then, we use 
those values of coupling constants to compare the magnetiza- 
tion curve of our model with the experimental one'"*. 

We use the triplet gap E^T (k = 0) = 26.7 K observed in 
inelastic neutron scatteringi^, the lower critical field H^i = 
10. 3T (or 15. OK if g = 2.17 is used), which marks the on- 
set of magnetization, and the saturation field''* Hc2 = 30. IT 
(43. 7K) as the experimental input. 

The triplet gap has been calculated in sec.|lll]and are given 
by eq.( |D5l ) or iD6i . In sec. [V] we have obtained the critical 
field iJci (eq.(i42b) and 11^2 (eq.(i45l)). We compare these re- 
sults with the experimental ones to determine two exchange 
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couplings Jl, J2 and the distortion parameter A. The result is: 

Ji = -140K, J2 = 87K, A = 0.46, (46) 

where we have used (ID6b for the excitation energy of the 
triplet. The magnetization curve for the ratio Ji/ J2 ~ —1.6 
and the distortion A = 0.46 obtained above is shown in 
FIG. [TtI (see FIG. [TS]). This curve is similar to that obtained 
in the high-field magnetization measuremen^i^ except for the 
little convex down character 




5 10 IS 20 25 30 

FIG. 17: Magnetization curve obtained from i39i by using the pa- 
rameter set l l46t . 



However, a remark is in order here. Recent NMR 
experiments^S suggest the displacement patterns of Cl~ which 
yield different magnetic interactions from what have been 
assumed here. In particular, the system does not have 
any explicitly tetramerized structure (see F1G[T]), although 
(CuCl)LaNb207 has period 2 both in the a- and the b direc- 
tion. Therefore, our results should not be taken literally. In- 
stead, our plaquette model should be thought of as one of the 
simplest Hamiltonians realizing the BEC of magnon bound 
states which is applicable to much wider class of systems in- 
cluding our simple J1-J2 model. 



VII. CONCLUSION 

Motivated by the recent discovery of a new two- 
dimensional spin-gap compound (CuCl)LaNb207, we have 
studied spin-1/2 J1-J2 model with a plaquette structure. For 
the small inter-plaquette interactions, i.e., for small A, there 
exists a finite spin gap over the spin-singlet ground state. 

We have computed the excitation energy of the triplets and 
the quintet in section in two different methods. If the 
gap of the lowest excitation closes, the corresponding particle 
condenses and a phase transition occurs from a paramagnetic 
phase to magnetically ordered phases. For the case of ferro- 
magnetic Jl considered here, we have two possibilities. For 
relatively small IJ1I/J2, the triplet particles (p and q) con- 
dense and generically we may expect CAF appears after the 
condensation (see FlG.fTOb. 

For larger values of | J1I/J2, however, the quintet excita- 
tion matters and we may have various phases. In the situation 
of relevance, we have either a usual ferromagnetic phase or 
a less conventional spin-nematic phase. One of these phases 



is selected by magnetic interaction among the quintet parti- 
cles. We have derived an effective Hamiltonian governing 
the magnetic part by using the second-order perturbation and 
mapped out the magnetic phase diagram (see FIG. [T2] |. A 
mean-field calculation predicted a finite window of the spin- 
nematic phase (green region in FlG.fTSli in agreement with re- 
cent numerical result&i!-' obtained for A = 1. From the proper- 
ties of the condensing particle, we found the nematic order for 
—2.33 < J1/J2 < —1.91 in the homogeneous (A = 1) J1-J2 
model in section HV] We remark that our effective Hamilto- 
nian is closely related to that for F = 2 cold atoms in optical 
Iattices213(i3i. 

We have studied the magnetization process in section[V] In 
the region of interest, magnetization is carried by spin-2 parti- 
cles, which should be identified with a tightly-bound magnon 
pair {magnon molecule),and we have constructed an effective 
hardcore boson (or, pseudo spin-1/2) model for these spin-2 
particles. By using a mean-field ansatz, we have determined 
the ground state of the above effective Hamiltonian as a func- 
tion of the external field h. We have found three different 
phases: (i) the normal BEC phase, (ii) the 'striped' supersolid 
phase and (iii) the 'striped' 1/2-plateau. In the normal BEC 
phase, the transverse magnetization (S^) vanishes unlike the 
conventional BEC in the spin-dimer model^. 

We have compared the results obtained for our Ji- J2 model 
with the experimental data of (CuCl)LaNb207 in section 
[VTl Although we have found that our model could qualita- 
tively explain the magnon gap in the inelastic neutron scat- 
tering experiments '^^ and the magnetization curveJi, the struc- 
ture suggested by NMR measurements^ is inconsistent with 
our tetramerized Ji- J2 model and this agreement should not 
be taken literally. Nevertheless, we hope that our scenario 
'molecular spin-BEC based on a simple J1-J2 model will 
capture the basic physics which underlies the magnetism of 
the compound (CuCl)LaNb207. 
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APPENDIX A: DIAGONALIZATION OF HAMILTONIAN BY 
BOGOLIUBOV TRANSFORMATION 

For convenience, we briefly summarize the method of Bo- 
goliubov transformation. We want to diagonalize 

H = Y.'^lA(k)vi., (Al) 

k 



14 



where 



vk = I Pk , gk , p^k ' 9-k 



(A2a) 



/ a(k) fe(k) c(k) d(k) \ 

-(^) ^ J!S =4'5 'SS ^ <-^' 

V&(k) rf(k) &(k) e(k)/ 

Now p, q are boson operators, and A(k) — A{~'k). We intro- 
duce Bogoliubov transformation 



ikVk = vj,, vj, = ( , , , q'\ 



where L is 4 x 4 real matrix, and 



(A3) 



Ll.ik) = L33(-k), LUk) = L34(-k) , 
Lt3(k)=L3i(-k), LtJk) = L32(-k) , (A4) 



[p'l 9'^] = -^11-^21 + L12L22 — L13L23 — ^14^24 = , 

\p'\p']=Ll,+Ll,-Ll,-Ll, = -l, 

(A5) 

where the summation over repeated indices is implied. By 
using 



and 



10 
10 
0-10 
-1 



e = 1 fori = 1,2 
e ^ -1 fori = 3,4 



(A6) 



(A7) 



the condition iA5\ can be rewritten as 



(A8) 



where 1 is a 4-dimensional column vector and the summation 
over i is not taken in (IA8b . g^i, can be considered as a metric. 
The condition (IA4b can be rewritten as 



r(r,2)(k) 



1(3,4) (-k) = 



(A9) 



where 1.; = ( r""^ , r^"^ 



and r is a 2-dimensional col- 



umn vector We denote A in (IA2bl i as and A = gA 



A^^ = g'^^Api,. Then, with regard to eigenvectors satisfying 

A^^M^ = aui^ , A^lu'j^ = bu'^, we obtain 



= T'^a^iA^^u'(^ , 



(AlO) 



since is symmetric. Therefore, if a 7^ 5, Ua^u'^ = 0, 
i.e. eigenvectors of a different eigenvalue are orthogonal each 
other. We define 



where = 1 ~ 4) are eigenvectors of A^]^, and 

1/2 _ 1 1/2 _ 1 
^(1,2) — -L' i(3,4) — ■ 

Then, L'gL'"^ = g. Therefore, 

L'^gL'g = I , 
where / is identity matrix. Now we can write 



A 



Ai A2 
A2 Ai 



(All) 



(A12) 



(A13) 



(A14) 



where Ai^2 are 2 x 2 matrix, and A(k) = A{—'k). Therefore, 
we can take = 1 ~ 4) satisfying iA9i . Defining fti as the 
eigenvalue of = 1 ~ 4), this leads to 



fli — — r^s, ^^2 
From ( |A13t , ( lAll i reduces to 



(A15) 



= V^L'gAL'^V 

= v'^L'.g ( fiil'i VL2V2 ^^313 ^^4l4)v' 

















^^2 



















V 











(A16) 



where v' = gL'gv, A — gA, and k is omitted. Comparing to 
( |A3t , we obtain 



L = gL'g . (A17) 
and ( |A9t . The eigenvalues of A (see eq. 



This L satisfies 
(IA2bl )) are given by 



r2(±.±) = ±[a — C — ff' 



'2 ' 



±{{-a^ + ^ + S-e^Y 



+ 4(a - d){c - e) (-46^ + ac + cd + ae + de)}^- 



(A18) 



We note that 

r^i 2 are the eigenvalues of the eigenvectors whose norm is 
positive, 

4 are the eigenvalues of the eigenvectors whose norm is 
negative. 
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and ( IA15I ) is satisfied. Moreover, using the boson's commuta- 
tion relation, H reduces to 

H = Y,{^i{^)pIA + ^2{\^)4^\.} + Eg , (A19a) 



{wi(k) +tj2(k)} 



(A19b) 



where uji,2 — "2.^1.2- 



APPENDIX B: HOW TO PARAMETRIZE GENERAL SPIN S 
STATES 

In this section, we briefly summarize the method of 
parametrizing arbitrary spin-2 states used in a mean-field cal- 
culation of section HVl The method is based on a geometrical 
representation of the spin-S" states used by Bacry-- and Bar- 
nett et al:^. Since our model has rotational symmetry, the 
mean-field energy has a trivial degeneracy with respect to the 
global rotation of the spin states. To mod out this degener- 
acy and find only essentially different solutions, this geomet- 
ric method is quite efficient. 

First we introduce the maximally polarized spin-S* state 
{spin coherent state) which is pointing the direction of 



(cos sm sm (p sm f/, cos ( 



i.e. (S-f2)|r2) = S\d). If we introduce the Schwinger bo- 
son operators a+ (a_) which destroys a spin parallel (anti- 
parallel) to the z-direction, the operator which creates a spin 
parallel to the r2-direction is given by 



(Bl) 



where 



■ X ■ ■ X ■ 4> 

u — e~'^e~' 2 cos - , u = e~'^e' 2 sin - (B2) 
2 2 

and X is an arbitrary gauge function. By using wt, the coherent 
state jO) can be written simply as 

1 



:(^t)2S|0) 



2S 



(B3) 



p=0 



where the combinatorial symbol 2sCp is defined by 2sCp = 

{2sy./ {{2s-py.pl). 

Next, we introduce a complex number C = {u/v)* ~ 
e"^ cot I and the corresponding unnormalized ket \C): 



2S 



p=0 



(B4) 



We note that the vector O rotates on the unit sphere S^, when 
the SU(2) rotation operator D acts on We denote an ar- 



bitrary spin-5 state by \A) ~ X]p=o^pl'^ ■ P ~ S). Then 



it is convenient to introduce the following 'wave function' 
which is in a one-to-one (except for an unphysical overall 
phase factor) correspondence with \A) under the condition 

2S 



p=0 



(B5) 



2S 



A2sl[{C-a,) 



where a^s are the 25 roots of (C) = and are parametrized 
as ai = e**^* cot y- If the degree deg of the above polynomial 
is smaller than 25, (25— deg) roots of P5 {() are at the infinity 
[di = or the north pole). Since the stereographic projection 
uniquely maps a set of 25 complex roots {ai} onto a set of 25 
points on a two-dimensional sphere S^, we can parametrize 
arbitrary spin-5 states by specifying 25 points on a sphere. 

If A2S = 0, the limit A2S 0, aj = 0{1/A2s) for any j 
must be taken {6j 0). In the case of spin-2, A'^ = Ai/A2s 
is given in terms of four complex numbers {oii] by 

A'q = aia2a3ai , 



A = - 



A' 
A' 



V6 



ai + a2 + as + a^ 



(B6) 



and hence the coefficients {Ai} read 



A,; 



-.A[ (z = 0,...,4), (B7) 



where (f) is the phase of A4. Therefore, as has been described 
above, arbitrary spin-2 states are parametrized by a set of four 
unit vectors and an overall phase factor. The rotational sym- 
metry enables us to further reduce the number of free param- 
eters by fixing ai and a2 as: 



ai = 1, aa = e"^\ a^. = e"^' cot y , 04 = e^^ cot -j- . 

(B8) 

Equations (|B6l l- (|B8l ) express arbitrary (except for global rota- 
tion) spin-2 states in terms of five free parameters. 



APPENDIX C: RELATION AMONG THE EXPECTATION 
VALUES OF SPIN-S OPERATORS 

There exists a simple relation among the expectation val- 
ues of spin-5 operators. By spin-5 operators, here we mean 
all independent (traceless) polynomials made up of the usual 
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spin-S* operators S. The spin 1 case has been considered by Using this relation, we obtain 
Chen and Levy^ in the context of spin-nematic order An 
arbitrary spin-S ket is written as 



Z2 



\ Z2S+1 / 



(CD 



^J—(4z,)i4,A (C4) 



where \zi\^ = 1. It is convenient to consider the Lie 

group SU(2S'+ 1) which naturally acts on the above {2S+1)- 
dimensional space. Let us denote the generators (a = 
1, . . . , (25 + 1)2-1) of SU{2S + 1) and normalize them as 



25+1 ' 



where the summation over repeated indices is implied. In 
spin-1/2 {SU{2)) case, can be written as ia", where a 
is Pauli matrix. Therefore, this relation can be written as 



(C2) 



Then, they satisfy 



{i,j,k,l = 1, 



,25+1) 



(C3) 



J 



APPENDIX D: EXPRESSIONS OF OMITTED EQUATIONS 



Section HITB] 

The elements of the second-order hopping matrix (eq. ( |24] |) are given by: 



a(k) 



-A|/+(k)-A2|/+(k) + A2 



-4Jf + SJfJa + 24Jf J| - 25JfJi + 28Ji J| - 28J| 

8J2(Jl-2J2)(Ji- J2)(Jl+2J2) 



J2 

4 



6(k)^-A^/_(k)-A2:^/_(k) + A2 



2 



-4Jf + 3JfJ2 + 24 Jfj| - 25JfJi + 28Ji J| - 28 J| J2 

T 



8J2(Jl-2J2)(Ji- J2)(Jl+2J2 



/± (k) = cos kx + cos ky + cos(fc^ ± fcy ) . 
The second-order energy shift for the singlet ground state is calculated as(see eq. (IZSll): 

3A2 (2 J2 + 3 J|) 



AS, = 



8J2 



The excitation gap of triplets from the second-order perturbation is given by: 
At 



^ ^3 2Jf-3J2J^ + J|jf„2j|j? 244^1-24^ 

*^ ^ 2' 8J2(J2- Jl)(2J2- Jl)(Jl+2J2) 



The modified excitation gap of triplets which is free from the pole Ji = 2 J2 is given by: 

At,„,od = E~^^^ik = 0) = J2 - A^J2 + A2 



3 2 25Ji- 65 J| 
2'^^ + ^ 144J2 



The excitation energy of quintet from the second-order perturbation is given by: 

'' -14Jf + 4J2 Ji + 15J|Ji - 6 J| {Jf + J|) 



£;g(k) = J1 + 2J2 + A-' 



8 Ji (2 Ji - J2 



4Ji 



(cos fcj. + cos ky ) 



Section HV] 

The parameters of the effective HamUtonian ( |32] | where the quintet condenses are given by: 

1 , _ _ A^ (59Jt + 78J2 Jf + 60J| + 26Jf Ji + J|) 



>/qui — ( Ji + J2) — 
o 



576Ji (3Ji2 +4J2J1 + J|) 



(C5) 



(Dl) 

(D2) 
(D3) 

(D4) 
(D5) 
(D6) 
(D7) 



(D8) 
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J, 



\,h AV| (l33Jf + 76J2 Ji + 7J| 



qu2 



16 



2304 Ji (3 J/ + 4J2 Ji + J|) 



(D9) 



qui 



2^^Jt + 174 Ja^^ + 120 J|Ji2 + 10J|Ji - J| 
2304 Ji {iJl + 4J2 Ji + J|) 



K„u_2 = — A' 



11J| + 68J1 J| + lilJlJl 
9216 Ji (37? +4J2J1 + J|) 



(DIO) 



r(l) _ 

-'^qul ~ 

r(4) _ 

''^qul ^ 



576 



>^'^{Jl + J2) A2) _„ .(3) _ 
-'^qul ^ 'J' ^qul ^ 



AV| ,(6) 
4608 Ji ' 'I"! 



AV, 



1152 



2 ^(5) 
1 ''^qul 



1152 ' 

AV|(J2 -3Ji) 
4608Ji ( Ji + J2) 



(DID 



L% = -2l';^i, forz = l,...,6. 

SectionlVl 

The parameters of the effective Hamiltonian iJeff ( l36b in the magnetization process are given in a series in A by: 



(D12) 



- M 
Jcsi = A 



A2 (-14Jf + 4J2 + 15 J| Ji - 6J|) 



8J1 (2Ji - J2) 
Ji + ,2 i^Jf + 11-^2 J? + 2J| Ji^ - 13J| Ji + 4J|) 



+ A^ 



16 Ji (2 Ji - J2) J2 



, -^J2 , .2^2(6Ji'-13J2Ji+4j; 
4 32Ji (2Ji - J2) 



(D13a) 
(D13b) 
(D13c) 



On top of them, we have three-body (or, three -plaquette) interactions: 



r(l) _ y2 (-^l +-^2) r(2)_n .(3) 



(4) 



(^1+^2), i^ff 



{J1 + J2), ifff' =0, i^ff =^ 



(6) 



\2-'2 



(D13d) 



The parameters necessary for the mean-field energy ( |39] l in the external magnetic field depend on the sublattice structures 
assumed in the calculation and are given as follows. 



1. In the case of "checkerboard" sublattice: 



4 - Jcffl + Jem + ^^eff + ^cS + ^off + 4^off ' 



(6) 



(1) 



"1 = n ( JcSi - L^^fl ) , Q;2 = - ( 2 Joff2 ~t i^^g 



-(1) 



4lS - 4L 



^cff 



(4) .(6) 

off ^off ; ' 



L 



(6) 
off ■ 



(D14) 



2. In the case of "striped" sublattice: 



£, = JcHl + JcS2 + ^^off 



(1) , r(3) , r(4) 



(6) 



-'^cft + ^cff + 4-'^oft 1 



ai = ^ ( Jcffl + 2 Jcff2 - L^ff - 4ieff - 4icff - 2ieff j , «2 = ^ i Jcffl + i^ff ; , 



(4) 



(6) 



(3i /32 = 2, 71 = ^ (iiff' + + 8L;- + 2L-^ j , 72 = ^^eff 



(3) 



(4) 



(6) 



(1) 



(6) 



(D15) 
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